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Abstract 

The equation for temporary evolution of spectral energy of collisionless Alfven waves is derived in 
framework of weak turbulence theory. The main nonlinear processes for such conditions are induced 
- scattering and two quantum absorption of Alfven waves by thermal ions. The equation for velocity 

t*^ ■ distribution of thermal particles is derived that describes diffusion in momentum space due to this 

. nonlinear processes. Comparison is done with the results of another authors. Results obtained are 

qualitatively differ from the ones obtained for the case of Alfven waves propagation along mean 

magnetic field. 
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Introduction 

It is well known that Alfven waves are very important agent of cosmic plasma. It is really 
believed now that main part of energy of interplanetary and interstellar magnetohydrodynamic 
\ (MHD) turbulence is contained in Alfven waves. This is because other types of low frequency 

^ ■ waves - ionsound and magnetosonic waves are strongly damped due to linear Landau damping. 

■ For Alfven waves such a damping is relatively small and nonlinear effects should be taken 
O , into account in order to investigate dissipation of MHD turbulence. It is well known now 
Q \ that the main nonlinear processes for such waves are induced scattering and two quantum 
absorption by thermal ions. Corresponding damping rates were calculated by many authors 
Qh' y^'"^]- Nevertheless totally comprehensive study is absent now. Firstly induced scattering of 

Q ■ Alfven waves was considered in paper [1] but for /5 -C 1. Here (3 is the square of the ratio of 
^ • thermal ion velocity vn = ^Ti/rui (Tj - is ion temperature in energetic units, mj is ion mass) 
p3 ■ and Alfven velocity Va = -Bo/V^vrpo (-Bo is the mean magnetic field strength, po is the mean 
plasma density). In papers [2-5] the case of wave propagation along mean magnetic field was 
^ ! considered. At last in papers [6,7] the oblique propagation of Alfven waves was considered but 
^ I random magnetic field component along mean magnetic field arising from magnetic pressure 
of Alfven waves was not taking into account. As will be shown in this paper the last effect 
strongly changes resultant damping rates. In particular, the interaction of waves with the same 
signs of wavevector components along mean magnetic field is impossible. 
Basic equations. 

Kinetic approach is the most adequate method in nonlinear plasma theory. In framework of 
weak turbulence it assumes the expansion of velocity distribution of thermal particles in series 
on powers of random electric and magnetic fields and calculation of nonlinear electric currents 
which are substituted in Maxwell equations. Such a procedure applied for Alfven waves results 
in cumbersome expressions containing multiple series of Bessel functions (cf. [1]). If the case 
of magnetized thermal particles is considered this functions should be expanded on powers of 
theirs small argument. These calculations are very tedious. We shall use another approach here 
that seems also more attractive from physical point of view. If wavelengths are much larger 
then thermal particles gyroradii it is possible to use equations averaged on gyrorotation. This 
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means that we shall neglect by all dispersive effects related with finite particle gyroradius. The 
system of equations can be written as follows [8] : 
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where h = B/B is a unit vector along magnetic field B, U£;=c[E x B]/5^ is electric drift 
velocity, i^y is electric field component along magnetic field, p and u are plasma density and 
velocity, Fa{v\\,v±), qa, ma are velocity distribution, electric charge and mass of thermal species 
a. Perpendicular and parallel pressures in Eq. (1) are given by the formulae 
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Guiding center equation (2) takes into account kinetic effects of thermal particles. For 
the case of frozen magnetic field that is considered here electric drift velocity is simply 
perpendicular to magnetic field component of medium velocity u. 
Expansion on powers of medium velocity. 

Let mean magnetic field Bq is in 2;-direction. Let us write all quantities in a form B = Bq + 5B, 
P = Po + ^P-i etc. and expand them in Fourier series: (5B = X^Bfeexp(ikr — iujt), Sp — 

k 

5pA: exp(ikr — icut) etc., k = (k,a;). All these Fourier components should be expanded on 

k 

powers of medium velocity due to presence of linear polarized Alfven waves. As for linear 
theory of Alfven waves Spk — 0, E\\k — 0, B^k — 0, u^k — 0, these quantities should be 
considered as second order ones. Performing expansion in Eq. (2) one can obtain second order 
SFak- Substitution of this quantity into quasineutrality equation (6) permits to find E\\k. As a 
result 
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For simplicity of notation we shall omit mean value of the product UkfUk" that should be 
extracted from correspondent product. Quantities (7Q{k',k") and in expressions (8) and 

(9) can be expressed in terms of velocity distribution of thermal particles: 
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Let us calculate 2;-component of random magnetic field B^k now. In order to do this one 
should expand Fourier transform of Eq. (1,3) to the second order and multiply Eq. (1) on kj_. 
The final result can be written as follows 
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Here quantity (72 (A;', k") is given by the expression 
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A{k) = klDo{k) + {clkl-uj')/vl 

corresponds to dispersion relation of magnetosonic wave A{k) = 0. 
The velocity Ca = 
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; — (P||o — P±o) I Po substitutes Alfvcn velocity Va in plasma with anisotropic 
pressure. We consider here the case of firehose stable plasma when the expression under square 
root is positive. The quantity Do{k) is given by the expression 
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Calculation of damping rates. 

Standard method of calculation of damping rates uses the expansion of equations (1-3) up to 
the third order of medium velocity. We use here more simple method that gives the same result. 
We shall derive equations for thermal particles and find probability of scattering. Practically 
one should derive quasilinear equation from Eq. (2). The only difference is that derivatives on 
velocity in Eq. (2) are multiplied on factors of second or higher order. This equation can be 
written as follows: 
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Substituting expression (12) for B^k and performing assemble averaging one can obtain 
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where l^(k) is spectral energy density of Alfven waves with dispersion relation c<;(k) =Ca|k^ 
It's normalized such that 

k 



First term in square brackets of expression (17) corresponds to two quantum absorption, the 
second one corresponds to induced scattering. The probability of two quantum absorption 
w(k', k", a;(k'),a;(k")) is given by the formula 
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Here ^p(k.'j_, k" ) is the angle between k'j^ and k^. The spectral energy density of Alfven waves 
evolves according to the equation 
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where damping rate can be expressed in terms of velocity distribution of thermal particles and 
probabilities of scattering and two quantum absorption [9,10]: 
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It is easy to see from expression (19) that waves with the same signs of kz can not interact. It 
is not so if one tend to zero k'j_ and k'^. This means that one-dimensional case is particular one 
(see Discussion). 

It is possible to transform damping rate (21) containing probability (19) to the more convenient 
for applications form: 
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and quantities Di{k) and D2{k) can be expressed in terms of velocity distribution of thermal 
particles: 
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It should be noted that damping rate (22) can be derived directly from expansion of Fourier 
transforms of Eq. (1-3) up to the third order of medium velocity but takes more tedious algebra. 
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Damping rates for Maxwellian plasma. 

Strictly speaking collision term should be added to the right-hand side of equation (17). It 
tends to make the velocity distribution MaxweUian. It will be so if frequency of coUisions 
is high enough. On the other hand frequency of collisions can be small enough for validity 
of coUisionless approximation. For Maxwellian plasma nonlinear damping rate (22) can be 
expressed in terms of function J+{x) [1]: 
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It contains principal value integral. 

For the sake of simplicity we shall consider pure hydrogen plasma. Damping rate is mainly 
determined by thermal ions. Input of thermal electrons is ^jmijme times smaller. The de- 
pendence of quantities DQ{k),Di{k), D2{k) on frequency and wavenumber is then reduced to 
dependence on x — 
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It was assumed here that ions and electrons can be at different temperatures Tj and Tg. We 
shall consider two extreme cases below. 
/5 < 1. 

For this case damping rate can be written in the form 
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where 6{x) is step function. 

The first term in expression (30) describes induced scattering of Alfvcn waves which can be 
interpreted as generation of ionsound wave by pair of Alfven waves that is absorbed by thermal 
particles. For the case considered scattering is differential, that is interaction of Alfven waves 
with close absolute values of kz is possible only. The second term describes two quantum 
absorption of Alfven waves which can be interpreted as generation of magnetosonic wave that 
is absorbed by thermal particles also. Expansion of spectral energy density of Alfven waves in 
the vicinity of —kz and change of summation to integration results in 
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First term in expression (31) corresponds to result obtained by Livshits and Tsytovich [1]. 
The only difference is the value of factor ol\ which depends on Tj/Tg in [1]. It is because in [1] 
simplifying assumption a; = in the denominator of integral (32) was used. The result obtained 
here is more correct and can be derived in framework of standard magnetohydrodynamics. 
/5> 1. 

For this case expression (22) can be reduced to 
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The damping rate (33) is determined by two quantum absorption that is again differential on 
absolute value of kz- Performing expansion similarly to previous case one can obtain 
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where 0^2 is given by the expression: 
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Numeric integration shows that this quantity practically does not depend on Ti/T^. Numeric 

value a2{Ti = Te) = 2.25. 

Discussion. 

The main result of this paper is the expression for nonlinear damping rate (22). Results obtained 
differ from results [6,7] where oblique Alfven waves were also considered. We take into account 
components of random field, corresponding to magnetosonic waves. As magnetosonic waves 
are strongly damped (especially in > 1 plasma, exceptions are quasiparallel and quasiper- 
pendicular propagation), these components should be considered as second order quantities 
and expressed in terms of medium velocity in Alfven waves. Simply speaking we treat these 
quantities similarly to parallel electric field component E\^k. This approach is correct if energy 
density of Alfven waves is not concentrated in wavevectors at small angles: 



< SB^ > /Bl < kl/kl 



(36) 
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For the opposite relation one-dimensional results [2-5] and results [6,7] are valid. The results 
of this papers can be readily obtained if one let B^k = in Eq. (16) or tend k±, k'j_ to zero in 
expressions (19) and (22): 



r(k) = L(k) Jm E cos^ ^(k^, k'^ 
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In particular interaction of waves with the same signs of kz is possible in this case. It is 
described by the second term of damping rate (37). Derivation of damping rate (37) shows 
that in order to obtain validity condition for approach used in this paper one should compare 
terms {kz ± k'^)'^ — (a; dz lu'^/cI and (k ± k')^, in denominator A(/c ± k') of damping rate (22). 
For the waves with the same singes of kz the fist term is equal to zero if dispersion relation 
for Alfven waves is used. It does not equal to zero if one takes into account nonlinear shift 
of frequency that is of the same order as nonlinear damping rate. This procedure results in 
condition (36). 

For /3 <S 1 linear damping of magnetosonic waves is smaller than for the case /3 > 1 and is due 
to Landau damping on thermal electrons [11]: 

7 = I^vt (38) 
V 8 ^ \kz\ rrii 

and approach used in this paper is valid for stronger condition: 

< 6B'^ > IrrieTe S 



«\H¥-'^Ti (39) 
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For the opposite relation linear damping of magnetosonic waves can be disregarded and results 
[1] are valid. Three wave interactions of Alfven and magnetosonic waves [11] should be also 
taken into account in this case. 
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